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Abstract 

In the functional integral approach to quantum gravity, the quantum configura- 
tions are usually treated to order U through a stationary phase approximation around 
the saddle point of the action where spacetime is flat. We show that from this point 
a "level line" in functional space departs, which comprises a family of static non-flat 
metrics with zero scalar curvature, depending on a continuous mass parameter. Fur- 
thermore, each of these metrics can be perturbed by an arbitrary function in such a 
way to still satisfy the condition / ^fgRd^x = 0. We thus find a set of zero- modes 
of the gravitational action which has non-vanishing measure in the functional space. 
These metrics will contribute to the functional integral as vacuum fluctuations, on 
the same footing as those near the saddle point. 

04.20.-q Classical general relativity. 
04.60.-m Quantum gravity. 
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1 Introduction 

The functional-integral approach to quantum gravity was initiated in the 1970s and early 
1980s by the monumental works of Hawking, Gibbons, Coleman and others pQ. Quantum 
effects at first order in K were computed in the stationary phase approximation through 
the zeta-function regularization technique. Hawking also tackled the problem of the un- 
definiteness of the gravitational action and showed that the action is unbounded from 
below, due to conformal modes which can be integrated away. The non-conformal part of 
the action was assumed to be positive semi-definite and this conjecture was then proved 
in several special cases [2]. Later, the un-definiteness of the action and the continuation to 
imaginary time were studied by Greensite, Mazur and Mottola and others [3]. Wetterich 
|4J found that the effective action is always un-defined in sign. 
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In the 1990s, Hamber and Williams implemented numerically a discrete version of 
Euclidean quantum gravity on a Regge lattice and computed non-perturbative quantities 
like the effective coupling [5 J . The Regge lattice approach also gives an insight into the cos- 
mological constant problem, showing that a space with average null curvature can emerge 
dynamically from a full quantum theory which includes strong (gravitational-only) vacuum 
fluctuations. 

In this work we shall consider pure gravity without any sources or matter fields, the 
action being given by the Einstein-Hilbert term plus the York-Hawking- Gibbons boundary 
term. We are thus looking at the vacuum state of quantum gravity. Our main result 
concerns the existence of a class of static non-flat metrics with zero action, which should be 
taken into account when computing physical quantities in the functional-integral approach. 

In our previous work [6] we found explicit expressions for these "zero-modes" in the 
weak field approximation. They owe their existence to the non-positivity of the Einstein- 
Hilbert lagrangian JgR. This feature makes possible a cancellation between positive and 
negative contributions to the integral of y^-R, arising from distinct spacetime regions. 

Here we obtain a generalization of the zero-modes to the strong field case. The 
situation can be understood through a finite-dimensional analogy. Imagine a system with 
an action which is not positive-definite. The stationary point corresponding to the classical 
theory will be a saddle point. The quantum configurations to order h span a small area 
around the saddle point. We can possibly find a "level line" (Fig. [1]) where the lagrangian 
has the same value it takes at the saddle point. Starting from the classical configuration 
at the saddle point, the system can evolve along this level line. Such lines do not exist, of 
course, for systems with positive-definite action. 

One should correctly object that level lines are subsets of null measure and do not 
contribute to an integration over the full space. But here the analogy fails to account 
for the complexity of a functional space. We shall prove that our strong-field zero-modes 
constitute a full-dimensional subset of the functional space of smooth, non-singular metrics, 
because they can be perturbed on an infinite-dimensional functional basis of the form 
{sin(27ma;), cos(27rmx)}. The oscillating weight factor exp(iS/h) in the functional integral 
(where S includes the boundary term) is constant over this set. 

The next question is: how can we take these configurations into account for the 
computation of physical averages? This is a very difficult task, and a problem in functional 
analysis which goes well beyond the traditional stationary-phase approximation. From 
the physical point of view, however, it seems clear that such strong vacuum fluctuations 
should make feel their presence in some way. One can envisage, for instance, some kind of 
gravitational Casimir effect, where the variation of external constraints on the fluctuations 
results in observable vacuum forces. 

In Sect. [2] we use a "virtual source" method to search for zero-modes of the action. 
In the context of wormholes physics jH [9] this method is also called "reverse solution of 
the Einstein equations" ; it consists in finding a source which generates a metric with some 
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Figure 1: Example of action level lines near a saddle point in two dimensions. A level line 
exists, where the action has the same value it takes at the saddle point. 

desired features. In our case, the desired metric is not easy to define, because it must 
obey the non-linear integral condition / y/g~Rd 4 x = 0. In the weak-field approximation this 
translates, up to terms of order G 2 , to the condition / T^S^d^x = 0. 

While in wormholes physics it is necessary to minimize and justify in some way the 
recourse to un-physical matter, here we are free to introduce any source, provided it satisfies 
the Bianchi identities. For weak fields it is sufficient to choose Tu = 0, leaving only the 
Too component, and build a virtual source satisfying the "dipolar" condition / T QO d 4: x = 0. 
This can be made, for instance, by a positive-density core surrounded by a negative shell 

®- 

In the strong field case, however, the pressure components of T^ u are essential in 
order to satisfy the Bianchi identity; from numerical solutions of the Oppenheimer-Volkoff 
equations, it turns out that the pressure tends to diverge in the negative-density shell, even 
beyond the point where the mass-energy density has been cut off. This way of looking for 
quantum configurations, by introducing unusual sources in the classical equations, seems 
therefore to be ineffective. 

In Sect. 02 we take a different, more straightforward approach. We consider the most 
general static symmetric metric and try to impose directly the zero-action condition. We 
obtain a closed solution in the special case of constant goo . In this case, the lagrangian 
density JgR has a relatively simple expression in terms of g rr , and we can invert this as a 
first-order differential equation for g rr in terms of y/gR. Setting R = 0, we obtain a family 
of non-flat metrics with zero scalar curvature, depending on a continuous parameter (the 
"level line"). Replacing yfgR with an arbitrary periodic function with zero average, we 
obtain the more general zero-modes mentioned above. 
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2 Search for zero-action modes with the virtual 
source method 



2.1 Weak field 

To first order in G, the field h^ u (x) generated by a mass-energy distribution T^ u (x) is given 
by an integral of the field propagator P pupa (x, y) over the source: 



M*0 = J ^yP^x^T^iv) (1) 

In Feynman gauge the propagator is given, with our conventions on the metric sig- 
nature, by 

Consider a static source in which only the component T 00 is non-zero. To order G, 
such a source automatically satisfies the Bianchi identities. The field is 

/+oo r 
dy / d 3 yT m (y)P lu/00 (x,y) 

2G f+°° f T 00 (y) 

= -(2r^o - VrVn) dlA, J A ^ _ yo)2 _ (x _ y)2 + %£ 

r T 00 (v) 

= 2G(2^ ^o - Moo) J d3 yj^Z^\ ^ 

The contracted form of the Einstein equations allows us to compute the scalar cur- 
vature, for any solution, as the trace of the energy-momentum tensor: 

R(x) = S7rGg^(x)T^(x) (4) 



In our case, it is straightforward to check that y g(x)g 00 (x) = 1 + o(G 2 ), and so the 
Einstein-Hilbert action computed for the solution of the weak-field equation is 

Szero-mode = — ^ / < ^xTqq{x) + o(G 2 ) (5) 

Provided the integral of the mass-energy density vanishes, the action of the field 
configuration is of order G 2 , i.e. practically negligible in several cases. One can give several 
examples of "dipolar" sources composed of positive and negative masses, either point-like 
or in the form of concentric shells, etc. [Bj. Some possible applications were mentioned in 

m 
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2.2 Strong field 



After giving an example of zero-modes in perturbation theory, let us look for dipolar 
fluctuations with spherical symmetry using the same virtual source approach, but in strong 
field. Our starting point is the well-known expression of the Schwarzschild metric (from 
[TO] , adapted to our conventions) 

dr 2 = e 2 ^dt 2 , - r 2 (d6 2 + sin 2 Odtf) (6) 

1 — 2m{r)G/r v ' 

The function m(r) represents the mass contained in a 3-sphere of radius r: 

Anr' 2 p(r')dr' (7) 



m r 



o 



where p(r) is the mass density in the local inertial frame. 

The function 0(r) appearing in the goo component of the metric is the solution of the 
system of coupled differential equations 

d(p(r) [m(r) + A7cr 3 p(r)]G 
dr r[r — 2m(r)G] 

m = - lP(r)+p(r /m (9) 

where p is the pressure. The second equation is called Oppenheimer-Volkoff equation and 
amounts to impose the Bianchi identities on the source. At the external boundary r ext of 
the mass distribution, <p satisfies the condition <p{r ext ) = (1/2) ln(l — 2MG/r ext ), where M 
is the total mass, M = m(r ext ). 

These relations are often employed in relativistic astrophysical stellar models. In that 
context, it is also necessary to assign a positive p; this is not necessary in our case because 
the source can be un-physical. We try to assign p in such a way to obtain metrics for which 
the integral / d A x^fgR is zero, by placing a negative-mass shell around the positive core, 
as done in the weak-field case. 

As a first example, let us consider a core with constant positive mass density p , and 
let r s be the corresponding Schwarzschild radius 



8irp G 



(10) 



This is the radius at which a singularity would arise, if the mass distribution was not 
cut before that point. Let us take an "inversion radius" ri nv and an external radius r ext 
such that r inv < r s < r ext , and define a density function p(r) which is equal to po up to 
r inv , becomes negative at r inv and then zero from r ext on: 

p(r) = p e(r inv - r)6(r ext - r) (11) 
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where e and 6 are the usual step functions. We shall try to define r ext in such a way 
to satisfy the zero-mode condition / d 4 x^/gR = 0. Note that r inv is arbitrary, provided 
finv < r s . If r inv is close to r s , the metric tends to become singular and the volume of the 
configuration can grow without limit. 

In order to write the action density we use the contracted Einstein equations (jlj) and 
the components of the energy-momentum tensor [10J 



T 00 = p/goo (12) 

T rr = p/grr (13) 

T ee = p/r 2 (14) 

T H = p /( r 2 sm 2 6) (15) 



We define the adimensional quantities 



s 


= r/r s 


p 


= P/Po 


p 


= PPo 


m 


= 2Gm/r s 




Rr 2 


R 
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We easily obtain in adimensional form 



v/l — 'fh/s 

The equations for <fi and p are, in adimensional form 

rh(s)/s + 3s 2 p(s) 



b'(s) 



(16) 
(17) 
(18) 
(19) 

(20) 



V9= , ; R = p-Zp (2i) 



(22) 



2s — m(s) 

p'(s) = -{p(s)+p(s)}<j ) '(s) (23) 

Next we make our treatment more general. Assuming that the virtual density p(r) 
decreases, in absolute value, when r increases, one can obtain several different metrics. 
The simplest case is that of functions of the form \p(r)\ oc l/r k , with k > 0. (Lobo [11] 
also employs a source density of this kind, in the context of wormhole solutions, with k = 2 
and without sign inversion. See also [12].) Eq.s (fTUj) . (TTTT) . ff20l and (122"]) are respectively 
replaced by 




(24) 

P( r )=Po(^j e(n - r)9(r 2 - r) (25) 
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R 



Rrl 



3 — k 

m(s)/s + (3 — k)s 2 p(s) 
2s — m(s) 



(26) 
(27) 



In order to avoid the density divergence in r = 0, we re-define p(r) as a constant 
from r = to, say, r = r s /2. It follows that the Schwarzschild singularity does not occur 
exactly at r s , but slightly on the right, at the point r = ar s , where a is the solution of the 
algebraic equation 

-- (28) 



a - a 3 k + 



3-2 3 " 



With this modification of the density function and this definition of a, we finally have 
in adimensional form 



p(s) = 9 (a — s) 



1 



a ( .s - - ) *-* 



+ 6(s-a)(-s- k ) 



(29) 



mis 




9(a 



+8(s - a) -s 



3-k 



k 



3 • 2 3 ~ k 



+ (30) 



The coupled equations (1221) and (1231 can now be solved numerically for given values 
of k. It is already apparent from the form of (123|) that the pressure can in general diverge, 
because the derivative of p is proportional to p itself, when the density p is small. In the 
usual astrophysical models, p goes rapidly to zero when r increases, and then p vanishes, 
too, at the star surface; in this way, p stays zero outside the star, as it is physically obvious. 
One would therefore expect that R = p outside a certain radius, and that the zero-mode 
condition can be satisfied by cutting p. 

In the present case, however, the situation is more tricky. The numerical solutions 
show that the pressure actually decreases quickly in the positive density core, but keeps 
non-zero in the outer negative-density shell. If the mass density is cut at r ext , the pressure 
begins to grow exponentially in the negative direction with a strange self-amplification 
effect (Fig. [2]), and the analogy with real sources is completely lost. According to (12T1) . it 
follows that the scalar curvature diverges, too, and there is apparently no way to obtain a 
zero- mode through this method (Fig. [3]). We have tried several values of k, in the range 
< k < 2. 



3 Zero-action modes with constant goo 

Since the virtual-source method appears to be ineffective, we look now for zero-action 
configurations in a way which is mathematically more straightforward and does not employ 
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Figure 2: Solution of the adimensional Oppenheimer-Volkoff equation (|23|) with k = 0.1, 
a = 1.00234 (eq. (|28|) ). density inversion radius 0.999a, density cutting point s — 1.1, initial 
condition p(0) = 10. The pressure decreases quickly but does not vanish, and diverges after 
the density cutting point. 




Figure 3: Adimensional Lagrangian density ^fgR (eq. (12T1) ) with the same parameters as 
in Fig. [2]. The scalar curvature diverges due to the pressure component, and the zero-mode 
condition can clearly not be satisfied. 
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the classical concepts of sources and solutions of the Einstein equations. We consider the 
most general static spherically symmetric metric and try to impose directly the zero-action 
condition / d^x^fgR = 0. Following Weinberg [13J we write this metric as 

dr 2 = B(r)dt 2 - A{r)dr 2 - r 2 {d6 2 + sin 2 6d<f) 2 ) (31) 

where A and B are arbitrary smooth functions. 

We add the requirement that outside a certain radius r ext , A and B take the 
Schwarzschild form, namely 

, . / 2GM\ t / , ( 2GM\- 1 „ , . 

B(r) = fl — J ; A{r) = ( 1 — J for r > r ext (32) 

This requirement serves two purposes: 

(1) It allows to give a physical meaning to these configurations, seen from the outside, 
as mass-energy fluctuations of strength M. For r > r ext their scalar curvature is zero. 

(2) The York-Hawking- Gibbons boundary term of the action is known to take in this 
case the form / = — M J dt pQ. It is therefore constant for our perturbed solution with the 
same M and L ^ (see below). 

Weinberg gives an useful expression for the Ricci tensor as a function of A and B: 

B " ^ a 

rr = 2B~l {a+) ~~r 

Ree = -l + ^(-a + b) + 
R(f><f> = sin 2 6Rgg 

r H r(a + 6) - 

2A 4A K ' 

where a = A' /A, b = B'/B and all functions depend only on r. 

In order to obtain the curvature scalar R, we contract R^ v with the components of 
g^ u . (Note that Weinberg's convention is dr 2 = —g^dx^dx", therefore g rr = A, g tt = —B, 



Rtt 





(33) 


1 

A 


(34) 




(35) 


B' 


(36) 


rA 



r 2 , Qm = r 2 sin 2 9.) We obtain 



fl= _|£ + ^ + 2 ^ (37 ) 

Given any two smooth trial functions A and B, it is straightforward to compute WgR- 
The problem is, that we do not have any control on the outcome, because the expression for 
R is complicated and highly non-linear; so imposing the integral condition J d A x^/gR = 
appears to be exceedingly difficult. 

We do find a set of solutions, however, if we make the drastic simplification goo = 
B = const. The scalar curvature multiplied by the volume element becomes in this case 

'rA' 



y/gR=-tocy/\BA\ f_ + l-_j (3* 
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Apart from the Newton constant G, L is the lagrangian density of the Einstein-Hilbert 
action, computed for this particular metric. 

Let us fix arbitrarily a reference radius r ext , and introduce reduced coordinates s = 
r/r ext . Define an auxiliary function a = A^ 1 . Regarding L(s) as known, eq. (|38|) becomes 
an explicit first-order differential equation for a 



1 a L\ \a\ 
a' = ---+ V ' ' (39) 



s s 



57TSi 



The boundary conditions (1521) are written 

B(s>l) = (l--); A(s>l) =[!--) (40) 




where M is a free parameter, the total mass in reduced units: M = 2GM/r ext . In the 
following we shall take M < 0, in order to avoid singularities. For r < r ext , we have 
B = B(l) = 1 — M. 



It is interesting to note that putting L = in eq. (1391) we can easily find an exact 
solution, ie a non-trivial static metric with R — 0. Namely, if 1 — a > 0, then a = 
1— e const / s, which does not satisfy the boundary condition; if 1— a < 0, then a = i-\- e const / S; 
implying e cost = —M. The resulting g rr component has the same form on the left and on 
the right of s — 1, namely 

9rr =(l + (41) 

while g 00 is constant and equal to (1 + \M\) for s < 1, and is equal to (1 + \M\/s) for s > 1. 

For large \M\ and small L, the last term in eq. (|39|) is a small perturbation. Since a 
never diverges and a -1 does not appear in the equation, the perturbed solution is not very 
different from ( 14TI) (see below). For values of \M\ of order 1 or smaller, the equation can 
be integrated numerically. If we choose a function L(s) with null integral on the interval 
(0, 1), we obtain a metric which is a zero-mode of the action but not of the lagrangian, like 
those we unsuccessfully looked for through the virtual source method. Take, for instance, 
L(s) = Lq sin(27rns), with n integer (Fig. HJ). 

In conclusion, we have found a family of metrics with null scalar curvature, depending 
on a continuous parameter M, which in our finite-dimensional analogy represent the level 
line R = 0. Furthermore, we have built a set of metrics close to the latter, by solving eq. 
(|39l) with L arbitrary but having null integral in the interval (0, 1). These metrics do not 
have zero scalar curvature, but still have null action. They make up a full-dimensional 
subset of the functional space; in order to show this, let us re-write eq. f l38l) in operator 
form as follows (remember that a(s) = g^{ s ) an d B = g 00 = const.): 

L(s) = -^£-(sa'(s) - a(s) + 1) = Ta(s) (42) 
\a\ 
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Figure 4: The g rr metric components of: a zero-mode with M = —2 and null scalar 
curvature (lower line); a zero-mode with the same M and adimensional lagrangian L = 
8 sin(67rs) (upper line, with three peaks); a zero-mode with the same M and L = 4 sin(87rs) 
(intermediate line, with four peaks). 

Eq. (|39|) can be written in operator form as 

a(s) = T- l L(s) (43) 

Consider functions a(s) and L(s) in the interval (0,1). Let Eq be the set of the 
metrics a satisfying the boundary condition (|32|) . Let E\ be the image of E through T, 
namely T : E — > E\. Any L{s)eE\ can be expressed in Fourier series as follows: 

oo 

L(s) = [o n sin(27rns) + b n cos(2nns)] (44) 

n=0 

Eliminate from E\ the functions with a 7^ 0. The remaining ensemble, E 2 , contains only 
functions with null integral in (0, 1). We make the conjecture that 

(1) The ensemble E 2 is infinite-dimensional. 

(2) The back-image T~ l (E2) (contained in i?o) is infinite-dimensional. 

A rigorous proof of this conjecture is not trivial. The conjecture could only be false, 
however, if the operator T _1 was quite singular, projecting E 2 into a finite-dimensional 
space. But this in not the case, as can be seen by solving the equation explicitily in the 
perturbative approximation of small L. If M > 1, after replacing a and B in the square 
root with their unperturbed solution (T4"T1) . eq. (1551) simplifies to 

a' 

which can be solved to give 

a{3) = 1 J 



la L 
+ 



s s ans 



3/2 



(45) 



1 H K -^= ) ds + const. 



57T V S 



(46) 
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We see that in this case the operator T -1 is linear and its nucleus contains only the function 
L(s) = —8ttt/s. Our conjecture is therefore true. On the other hand, if M <C 1, eq. (|39|) 
takes the form 



Also in this case the operator T 1 is linear, its nucleus containing only the function L(s) = 
— ^7=, and our conjecture is true. 

4 Conclusions 

Although a generally accepted theory of quantum gravity does not yet exist, the aim of 
this work was to give some insight into its vacuum state. To this end, we considered the 
formal path integral over all metric configurations, weighed by the factor e lS ^ n . This factor 
oscillates in general very fast, and the most probable configurations are those which make 
the phase S/h constant or stationary. The metrics which make the action stationary are, 
of course, the classical configurations. A WKB-like expansion near these configurations in 
powers of h has been previously considered by Hawking and others in the context of the 
Euclidean functional integral approach. Usually, in quantum mechanics or quantum field 
theory the configurations with stationary action are local minima; so if the system leaves 
them, the phase of the weight factor starts to oscillate more and more quickly, until the 
one-loop approximation fails to be valid and a pure quantum regime is entered. In quantum 
gravity, on the contrary, the non-positivity of the action implies that it is possible to depart 
from the stationary point along configurations at which the action takes the same value 
it takes at the stationary point. These configurations actually have vanishing Lagrangian 
density. Each of them can be furthermore perturbed, or "deformed", in such a way as to 
obtain metrics which do not have zero Lagrangian at each point, but have zero total action 
due to the cancellation between positive and negative contributions from different regions 
of space. 

These zero-modes are very different from the classical solutions of Einstein equations 
and can not be obtained employing virtual sources or starting from wormhole-like metrics. 
It was actually to be expected that a full quantum feature of spacetime could not be 
recovered by "quantizing the classical solutions". The static and spherically symmetric 
zero-modes described in this work are probably just members of a much wider class. They 
are characterized by a large negative mass. They are not singular and outside a certain 
radius they look like Schwarzschild solutions. These metrics are physically realized as 
vacuum fluctuations and are therefore homogeneously distributed in space and time, but 
could possibly manifest their presence if some cause of inhomogeneity is introduced, as 




(47) 



whose solution is 





12 



it happens in electro magnetism for the Casimir effect. Finally, they could give a relevant 
contribution to the cosmological term. These latter issues will be addressed in future work. 
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